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Q_i| Abstract 

. An interesting discovery in the last two years in the field of mathematical physics 

has been the exceptional Xg Laguerre and Jacobi polynomials. Unlike the well-known 
classical orthogonal polynomials which start with constant terms, these new polyno- 
mials have the lowest degree I = 1, 2, . . ., and yet they form complete sets with respect 

t— I \ to some positive-definite measure. In this paper, we study one important aspect of 

these new polynomials, namely, the behaviors of their zeros as some parameters of the 

\q • Hamiltonians change. 

vo : 

in ■ 

H ! 1 Introduction 

o ■ 

T— I ■ 

The discovery of new types of orthogonal polynomials, called the exceptional Xg polynomials, 
has been the most interesting development in the area of exactly solvable models in quantum 
mechanics in the last two years [TJ |2j [3j HJ [5] . Unlike the classical orthogonal polynomials, 
these new polynomials have the remarkable properties that they still form complete sets 
with respect to some positive-definite measure, although they start with degree I polynomials 
instead of a constant. Four sets of infinite families of such polynomials, namely, the Laguerre 
type LI, L2, and the Jacobi type Jl, J2 Xg polynomials, with £ = 1,2,..., were constructed in 
[3] . These systems were derived by deforming the radial oscillator potential and the Darboux- 
Poschl- Teller (DPT) potential in terms of an eigenpolynomial of degree I (£ = 1,2,...). 
The lowest (£ = 1) examples, the XpLaguerre and Xx-Jacobi polynomials, are equivalent to 
those introduced in the pioneering work of Gomez-Ullate et al. [T] within the Sturm-Liouville 
theory. The results in [1] were reformulated in the framework of quantum mechanics and 
shape-invariant potentials by Quesne et al. [2j. By construction these new orthogonal 



polynomials satisfy a second order differential equation (the Schrodinger equation) without 
contradicting Bochner's theorem [6] , since they start at degree I > instead of the degree zero 
constant term. Generalization of exceptional orthogonal polynomials to discrete quantum 
mechanical systems was done in [7]. 

Later, equivalent but much simpler looking forms of the Laguerre- and Jacobi-type Xg 
polynomials than those originally presented in [3] were given in [3J. These nice forms were 
derived based on an analysis of the second order differential equations for the Xg polynomials 
within the framework of the Fuchsian differential equations in the entire complex x-plane. 
They allow us to study in-depth some important properties of the Xg polynomials, such 
as the actions of the forward and backward shift operators on the Xg polynomials, Gram- 
Schmidt orthonormalization for the algebraic construction of the Xg polynomials, Rodrigues 
formulas, and the generating functions of these new polynomials. 

Recently, these exceptional orthogonal polynomials were generated by means of the 
Darboux-Crum transformation [HE]. Physical models which may involve these new polyno- 
mials were considered in 1101. 



One important aspect related to these new polynomials, which was only briefly mentioned 
in |3] but has not been investigated in-depth so far, is the structure of their zeros. It is the 
purpose of this paper to look into this. Particularly, we investigate the behaviors of the zeros 
as the parameters of the polynomials change. 

The plan of this paper is as follows. In Sect. 2 we briefly review the forms of the 
exceptional polynomials. Sect. 3 and 4 study the behaviors of the extra and the ordinary 
zeros, respectively, of the exceptional polynomials as one of I and n increases while the other 
parameters being kept fixed. Sect. 5 presents analytical proofs that explain the movements 
of the extra zeros of the exceptional polynomials as n changes at fixed I. In Sect. 6 we 
consider behaviors of the zeros at large g and/or h. Sect. 7 summarizes the paper. 

2 Exceptional orthogonal polynomials 

Four sets of infinitely many exceptional orthogonal polynomials were derived in [3], among 
them two are deformations of the Laguerre polynomials, and the others are deformations of 
the Jacobi polynomials. A unified nice form of these polynomials was given in [I], in which 
these polynomials are expressed as a bilinear form of the original polynomials, the Laguerre 
or Jacobi polynomials and the deforming polynomials, depending on the set of parameters 
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A and their shifts d and a non-negative integer £, which is the degree of the deforming 
polynomials. The two sets of exceptional Laguerre polynomials (£ — 1, 2, . . ., n = 0, 1, 2, . . .) 
are: 

r fcfa; A + S)P n ( V ; g + i-1)- frfa X)d v P n (r ] ; g + i-1) : LI 
P^fajA)^ 0/ + J)6(i7; A + «)P„( W + ^ + , (1) 

[ +^;A)^P„(»7;</ + ^+l) :L2 

in which A = f g > and <5 = f 1 and 

9) d = f iSr^ to a = ( CS: 1 ? : 11 • ^ 

l ^ 2 fa) : L2 
The two sets of exceptional Jacobi polynomials {£ = 1, 2, . . ., n = 0, 1, 2, . . .) are : 

f(A + |)to A + «)P n fa; (/ + £ - 1, h + £ + 1) 
p , ^defl +(l + r / )^fa;A)^P n (r / ;^ + £-l,/ i + £+l) : Jl 

{ - (1 - 77)6(77; X)d v P n ( m g + £ + 1, h + £ - 1) : J2 

in which A = f (g, /i), (7 > 0, /i > 0, d = f (1, 1) and 

P e (77 , g > h > : Jl 

P £ 2 2 (77), h > g > : J2 

The new exceptional orthogonal polynomials can be viewed as deformations of the classical 
orthogonal polynomials by the parameter £, and the two polynomials £^(77; A) and 6(77; A + 8) 
played the role of the deforming polunomials. 

The zeros of orthogonal polynomials have always attracted the interest of researchers. In 
this paper we shall study the properties of the zeros of these new exceptional polynomials 
as some of their basic parameters change. 

In the case of Xi polynomial P^ n {j]\ A), it has n zeros in the (ordinary) domain where the 
weight function is defined, that is (0, 00) for the LI and L2 polynomials and (—1, 1) for the 
Jl and J2 polynomials. The behavior of these zeros, which we shall call the ordinary zeros, 
are the same as those of other ordinary orthogonal polynomials. We shall say more about 
these zeros in sect. SJ Besides these n zeros, there are extra £ zeros outside the ordinary 
domain. For convenience, we shall adopt the following notation for the zeros of the various 
polynomials involved: 

£f : zeros of £,(77; A + 8), k = 1, 2, . . . , £; (5) 



zeros of &fo;A), k — 1,2, . . . ,£; (6) 

7^f' n ^ : extra zeros of P^ n , k = 1, 2, . . . , £; (7) 

^.™) . Quinary zeros Q f p ln ^ j — 1,2, . . . ,n. (8) 

We emphasize that rf^ G (0, oo) for the LI and L2 Laguerre polynomials, and r/j 1 '^ G 
(— 1, 1) for the Jl, and J2 Jacobi polynomials. 

Figs. 1-10 depict the distribution of the zeros for some representative parameters of the 
systems, namely, n, £, g and h. From these figures one can deduce certain patterns of the 
distribution of the zeros as those parameters vary. We will discuss these behaviors below. 

3 Extra t zeros of Xi polynomials 

Here we discuss the location of the extra I zeros of the exceptional orthogonal polynomials, 
which lie in various different positions for the different types of polynomials. From our 
numerical analysis, we can summarize the trend as follows. 

The I extra zeros of LI polynomials are on the negative real line (— oo, 0). The L2 Xt. Q && 
polynomials have one real negative zero which lies to the left of the remaining |(£ — 1) pairs 
of complex conjugate roots. The L2 Xi :even polynomials have |£ pairs of complex conjugate 
roots. 

The situations for the Xi Jacobi polynomials are a bit more complicated. The Jl X^. ^ 
polynomials have one real negative root which lies to the left of the remaining — 1) pairs 
of complex conjugate roots with negative real parts. The Jl X£ :even polynomials have ~£ pairs 
of complex conjugate roots with negative real parts. The J2 Xi- ^ polynomials have one 
real positive root which lies to the right of the remaining — 1) pairs of complex conjugate 
roots with positive real parts. The J2 Xi- even polynomials have \t pairs of complex conjugate 
roots with positive real parts. 

One notes that the Jl and J2 polynomials are the mirror images of each other, in the 
sense r\ -H- —77 and g -H- h, as exemplified by the relation ^p{r];g,h) = (— 1)^^/ 1 (— 77; h, g) 
[31 0]. So the behaviors of the zeros of Jl Jacobi polynomials can be obtained from those 
of the J2 type accordingly. As such, for clarity of presentation, we shall only discuss the 
behaviors of the zeros of the J2 Jacobi polynomials in this paper. 
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Table 1: List of the zeros £r,£? and ^' n) for the LI Laguerre polynomials with g = 2, £ = 5 



and n = 0, 10, 20, . . . , 60 (k = 1, 2 . . . , £). It can be seen that when n = 0, r^' n " 0) = f^. As 
n increases, 77]^ ^ approach to £^ . 

^ : -22.4802 -15.2391 -10.1403 -6.2977 -3.3427" 

n = -22.4802 -15.2391 -10.1403 -6.2977 -3.3427 

10 -22.0686 -14.8767 -9.8314 -6.0505 -3.1698 

20 -21.8830 -14.7189 -9.7004 -5.9469 -3.0962 

f)f n) : 30 -21.7717 -14.6253 -9.6233 -5.8862 -3.0529 

40 -21.6954 -14.5617 -9.5711 -5.8452 -3.0237 

50 -21.6390 -14.5148 -9.5327 -5.8152 -3.0022 

60 -21.5951 -14.4784 -9.5030 -5.7919 -2.9856 

: -21.0456 -14.0274 -9.1375 -5.5071 -2.7824~ 



Table 2: Same as Table 1 for LI Laguerre polynomials with g = 8 and £ = 5. 



m . 

Sfe • 




-30.7592 


-22.3415 


-16.1499 


-11.2032 


-7.0462 




n = 


-30.7592 


-22.3415 


-16.1499 


-11.2032 


-7.0462 




10 


-30.4724 


-22.0859 


-15.9269 


-11.0165 


-6.9029 




20 


-30.3144 


-21.9474 


-15.8074 


-10.9169 


-6.8255 


-(e,n) 


30 


-30.2107 


-21.8574 


-15.7301 


-10.8525 


-6.7752 


40 


-30.1361 


-21.7928 


-15.6748 


-10.8065 


-6.7393 




50 


-30.0791 


-21.7436 


-15.6328 


-10.7715 


-6.7119 




60 


-30.0336 


-21.7046 


-15.5994 


-10.7438 


-6.6902 


M . 




-29.4106 


-21.1735 


-15.1488 


-10.3703 


-6.3968 



3.1 Behaviors as n increases at fixed i 

In all cases, we have 

p ifi (vA)^UvA + S). (9) 

This implies that the zeros of P^o coincide with those of ^(77; A + S), namely, ^ , k = 

1,2,..., e. 

At fixed i, all the fj^'^ move from £^ at n = 0, to ^ asn-f 00. This can be seen from 
Figs. [TH3]and in Tables HHZl We shall prove this result generally in Sect. |5j 
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Table 3: Same as Table 1 but for L2 Laguerre polynomials with g = 3 and £ = 4. 



At) 




-5.29007 ± 1.65310 i 


O ^ C\C\C\ O 1 V f\ V 1 OA 

-3.70993 ± 5.05130 i 




n = 


-5.29007 ± 1.65310 i 


-3.70993 ± 5.05130 i 




10 


-4.84198 ± 1.57129 i 


-3.25524 ± 4.78004 i 




20 


-4.71299 ± 1.54888 % 


-3.12839 ± 4.70776 % 


AW 
'Ik 


: 30 


-4.64523 ± 1.53732 % 


-3.06246 ± 4.67065 i 


40 


-4.60183 ± 1.52998 i 


-3.02046 ± 4.64713 i 




50 


-4.57100 ± 1.52479 % 


-2.99074 ± 4.63053 % 




60 


-4.54766 ± 1.52087 i 


-2.96828 ± 4.61801 i 


At) . 




-4.28361 ± 1.47684 % 


-2.71639 ± 4.47739 i 



Table 4: Same as Table 3 for L2 Laguerre polynomials with g = 10 and 

^ : -12.8111 -12.2115 ± 4.7185 \ -10.1329 ± 9.7965 \ 

n = -12.8111 -12.2115 ± 4.7185 i -10.1329 ± 9.7965 i 

10 -12.5476 -11.9465 ± 4.6639 % -9.8622 ± 9.6780 % 

20 -12.4210 -11.8198 ± 4.6384 % -9.7348 ± 9.6233 % 

r\f n) : 30 -12.3430 -11.7418 ± 4.6229 i -9.6570 ± 9.5901 i 

40 -12.2888 -11.6877 ± 4.6122 i -9.6032 ± 9.5673 i 

50 -12.2483 -11.6473 ± 4.6043 % -9.5632 ± 9.5504 % 

60 -12.2165 -11.6157 ± 4.5981 i -9.5319 ± 9.5372 i 

£jP : -11.8092 -11.2107 ± 4.5195 i -9.1347 ± 9.3702 i 



Table 5: Same as Table 1 but for J2 Jacobi polynomials with g — 3, h = 4 and £ = 4. 







1.56846 ± 2.10278 i 


3.00297 ± 0.91199 i 




n = 


1.56846 ± 2.10278 % 


3.00297 ± 0.91199 i 




10 


1.45201 ± 1.89890 % 


2.76834 ± 0.82626 i 




20 


1.42407 ± 1.85433 % 


2.71360 ± 0.80733 % 


'Ik 


30 


1.41139 ± 1.83435 i 


2.68882 ± 0.79884 i 




40 


1.40414 ± 1.82297 % 


2.67466 ± 0.79401 i 




50 


1.39944 ± 1.81561 % 


2.66550 ± 0.79088 % 




60 


1.39615 ± 1.81046 % 


2.65907 ± 0.78869 i 


S) . 




1.37745 ± 1.78118 i 


2.62255 ± 0.77624 i 
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Table 6: Same as Table 5 for J2 Jacobi polynomials with g = 3, h = 4 and £ = 5. 



m . 


1.19188 ± 1.85256 % 


2.38851 ± 1.21416 % 


2.83923 


ra = 


1.19188 ± 1.85256 % 


2.38851 ± 1.21416 % 


2.83923 


10 


1.11856 ± 1.68660 % 


2.22979 ± 1.11021 % 


2.64753 


20 


1.09936 ± 1.64851 i 


2.18998 ± 1.08600 % 


2.59983 


€ nY - 30 


1.09041 ± 1.63110 i 


2.17151 ± 1.07491 i 


2.57771 


40 


1.08522 ± 1.62106 i 


2.16081 ± 1.06852 i 


2.56490 


50 


1.08184 ± 1.61453 % 


2.15382 ± 1.06436 % 


2.55654 


60 


1.07945 ± 1.60993 % 


2.14890 ± 1.06143 % 


2.55066 


At) . 


1.06566 ± 1.58339 i 


2.12047 ± 1.04452 % 


2.51663 



Table 7: same as Table 5 for J2 Jacobi polynomials with g = 8, h = 9 and £ = 3. 







3.90615 ± 4.35051 % 


6.58770 


n 


= 


3.90615 ± 4.35051 % 


6.58770 




10 


3.74981 ± 4.16635 i 


6.32188 




20 


3.69527 ± 4.10323 % 


6.22948 


'Ik 


30 


3.66745 ± 4.07116 i 


6.18240 




40 


3.65057 ± 4.05174 % 


6.15383 




50 


3.63924 ± 4.03870 i 


6.13465 




60 


3.63110 ± 4.02935 % 


6.12088 



: 3.58151 ± 3.97238 i 6.03699 
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3.2 Behaviors as £ increases at fixed n 

The discussions in the last subsection show that fj^^ are sandwiched between £^ and ^ . 
Thus to know how fj^^ behave as £ increases at fixed n, we only need to study how the 
zeros and £^ flow as £ increases. 

3.2.1 LI Laguerre 

As £ changes to £ + 1, the zeros of £^(77; g + 1) and £^(77; g) decrease (move to the left), and a 
new set of zeros appear from the right. 

for jfe = 1,2,..., £- 1,£. 

We show these patterns for some representative parameters in Figs. H] and [5] 

3.2.2 L2 Laguerre 

For £ = 1, there is one real root each for ^(77; g> + 1) and ^(77; g), with £^ < £j < 0. 

For £ = 2, the above two roots bifurcate into two complex roots, with 9?<^ £ ) < 9R£^, > 

Generally, for even £, there are £ complex zeros with 

3tg><*ttP,\S!$\>\$*$\, k = l,2,...,£/2. (10) 

All are sandwiched between £^ and ^ . As an even £ changes to £ + 1 which is odd, all 
zeros move to the right with the real and the absolute value of the imaginary parts increased, 
and a new real zero appears to the left of all the complex zeros on the negative real axis. As 
£ increases further, the complex zeros move as described before, and the zero on the negative 
real axis bifurcates into two complex zeros, giving an even number of complex zeros. These 
patterns continue as £ increases. 

Figs. [6] and [7] show these behaviors for some selected parameters. For large £, these zeros 
distribute in a horse-shoe pattern. 
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3.2.3 J2 Jacobi 

For I — 1, there is one real root each for ^(77; # + 1) and £^(77; g), with £j > £^ > 1. 

For £ = 2, the above two roots bifurcate into two complex roots, with > 9ft£^, | > 

|^)|. 

Generally, for even £, there are I complex zeros with 

^l e) >^i E \mk ) \>m i k ) i k = i,2,...,e/2. (11) 

As I changes to i + 1 which is odd, all zeros move toward the y-axis, with the real parts 
decreased, and a new real zero appears to the right of all the complex zeros on the real 
x-axis. As I increases further, the complex zeros move as described before, and the zero 
on the real axis bifurcates into two complex zeros, giving an even number of complex zeros. 
The absolute value of the imaginary part of the complex zeros may increase initially, but 
eventually decrease as I increases. This pattern continues as I increases. 

Figs. M and [9] show these behaviors for some selected parameters. For large £, these zeros 
distribute in a horse-shoe pattern. 

4 Ordinary zeros of polynomials 

In the case of Xf polynomials Pi^ n {r\\ A), it has n zeros in the (ordinary) domain where the 
weight function is defined, that is (0, 00) for the LI and L2 polynomials and (—1, 1) for the 
Jl and J2 polynomials. The behavior of these zeros are the same as those of other ordinary 
orthogonal polynomials. 

4.1 Behaviors as n increases at fixed I 

This is guaranteed by the oscillation theorem of the one-dimensional quantum mechanics, 
since P^ n (r)\ A) are obtained as the polynomial part of the eigenf unctions of a shape invariant 
quantum mechanical problem. Explicitly, as n changes to n + 1, all zeros P^ n decrease, and 
a new zero appears from the right. Thus the n zeros of Pt^rSjl'-, A) and the n + 1 zeros of 
Pi,n+i(V) interlace with each other: each zero of Pt.rSjW A) is surrounded by two zeros of 
P^„+i(r?; A). 

Figs. show these behaviors for selected parameters. 
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4.2 Behaviors as £ increases at fixed n 

From Figs. 0]-[71 one sees that for LI and L2 Laguerre polynomials (whose zeros are positive 
in the ordinary domains), all the n zeros shift to the right as £ increases. 

For J2 Jacobi polynomials, the positive (negative) zeros shift left (right) as £ increases, 
i.e., they move toward the origin rj = 0. This is illustrated in Figs. |8]and[9j 

4.3 Additional observation for the LI case 

Using the well-known derivative relation 

d v Li a) (?l) = -Lt + i\v) (12) 

and 

L^ ) (v)-L^- 1 \v) = L i :lM, (13) 

we get 

Hence when n = £, the LI Laguerre is an even function of r], and its zeros are symmetric 
w.r.t rj = 0. 

5 Proof that f]^^ — > gj® as n — > oo 

As mentioned before, for n = 0, we have ff^'^ = £^ , as P n =o — 1- We shall show that as 
n — > oo, fil^ — > ^ . This amounts to showing that in this limit, d n P n dominates over P n 

5.1 LI and L2 cases 

We shall make use of the above derivative relation (II 2p and (Perron) Theorem 8.22.3 of [TT] . 
namely, 

L^\ V ) = -^-(-^-f-inf-le 2 ^^, aeR, rj e C\(0, oo), (15) 

2V7T 

which gives the asymptotic form of Ln(rj) for large n. For the LI and L2 cases, we have 
a = g + £ — 3/2 and g + £ + 1/2, respectively. 
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One finds 



n 



-Tj) 



(16) 



For large n with fixed rj, d v L^n\f]) dominates over L^(rj), and thus the zeros of Pe, n are 
determined by those of ^(rj; g) as n — > oo. 

5.2 J2 Jacobi 

For the asymptotic form of Pn a '^\r]) for large n, we shall make use of Theorem 8.21.7 of [TTj : 



/„2_i\-l/4 



v/2 



7m 



(17) 



ajel, r? G C\[-l,l] 



and 



d v P^\ V ) = l( n + a + P + 1)1^^%). 



One finds 



d n Pt 0) (v) {ri + a + (5 + l)\l n {V > (^T+T+v^F 7 !) 2 ' 



(19) 



Again, for large n with fixed rj, d v Pn (rj) dominates over Pn (rj), and thus the zeros 
of P^ n are determined by those of Ceiv'i 9) as n ~ * °°- 

6 Behaviors at large g and /or h 
6.1 LI Laguerre 



As g increases, we have 



ISA; I) l?fc I 



1% I 



I 

1% I 



(20) 



all increase. That is, all the zeros move away from the y-axis. This can be seen from Figs. H] 
and 

In fact for large g, we have ^(77; g + 1) ~ 6(77; <?)■ Hence 
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« Uv;gW^\v). (21) 

For 5» 3> 1, P^niV'i 9) approaches 

P^{m9)^Lf^{-r,)L^{ n ). (22) 

Thus the extra (^'™' ) ) and the ordinary (t]^' n ^) zeros of P^ n (rj]g) are given by the zeros of 
(—rj) and Ln (rj), respectively. 

6.2 L2 Laguerre 

As g increases, we have 9ft£^\ decreased, l^^l, |9?^ f | increased, and ^' n ' ) increased. 

This is easily seen from Figs. [6] and [71 That is, the zeros ££\ \ and hence fj^' n \ all are 
moving leftwards and away from the x-axis, while the ordinary zeros are moving towards 
the right. 

In fact for large g, we have ^(77; g + 1) ~ £,e{v'i 9)- Hence 



Using Eqs. (E.2), (E.10) and (E.9) of 0], we arrive at 



(23) 



Pt, n (mg) « eK^;^) + ^ + ^ + • (24) 

For g ^> 1, Pi n (rj;g) approaches 

P^jyJwL^^Lto+O^). (25) 

Thus the extra (fj^ ) and the ordinary (^ ,n ' ) ) zeros of P^ n (rj; g) are given by the zeros of 
h\ 9 ir\) and Ln +i \r]), respectively. 

6.3 J2 Jacobi 

As g, h increases, we have , ^Cfe , I^Cit |> I^Cfc I increased, as is evident from Figs. [8] 
and |9j The extra zeros fj^' n \ being in between these zeros, follow the same pattern. That is, 
the zeros ^\ ^ , and hence ff k , all are moving away from the x and y-axes. The ordinary 
zeros T)f will have their norm 1^ | decrease in general as g increases. Thus these zeros 
move towards the y-axis. 
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In fact for large g and h, we have (a = g + £ + ~, f5 = /i + £ — |) 

PiAv,9,h) « toP.fc) (^^^(^-(l-^^P^^) . (26) 
Using Eqs. (E.13) and (E.23) of [1], we arrive at 

^»(i75 0, *0 ~ <?, fc) [(a + n) P^ 1 ^ 1 )^) - £P^\v)] . (27) 

For g ^> 1 and h 1, Pe,n(v'i 9i h) approaches 

<?) « P^^P^^to). (28) 

Thus the extra (fj^ ) and the ordinary {r^' n ^) zeros of Pe jTl (j]; g, h) are given by the zeros of 
pi g e,h+i)^ anc [ p(s+t,h+£) ^ respectively. 

6.3.1 Additional observation: h^> g 

For l}> g, all zeros, i.e., £• , £• , ^f' n \ ?7fc £ ' ra ' ) , gather around 77 = 1. This can be understood 
as follows. From the series expansion of the Jacobi polynomials, Eq. (E.ll) of [I], 

T>(cP)f„\ - ( Q + X )" f 1 (- n )k( n + a + {3 + l) k ( l-V \ k ( Q \ 

" W " n! 2-j fc! (a + i) fe I 2 J ' (2Jj 

fc=0 

one sees that, for h ^> g, the absolute value of Pe, n {v', 9->h) is large near 77 = — 1 and small 
at i] = 1. Hence, in this limit, the zeros of Pi, n {j}\ g, h) distribute very near 77 = 1. We show 
this in Fig. [10] for certain parameters. 

7 Summary 

The discovery of new types of orthogonal polynomials, called the exceptional Xi Laguerre and 
Jacobi polynomials has aroused great interest in the last two years. Unlike the well-known 
classical orthogonal polynomials which start with constant terms, these new polynomials 
Pe, n (Vi ^) have the lowest degree £ = 1,2,..., and yet they form a complete set with respect 
to some positive-definite measure. Many essential properties have been studied in [I]. 

In this paper, we have considered the distributions of the zeros of these new polynomials 
as some parameters of the Hamiltonians change. The Xn polynomials Pg }n (r); A) has n zeros 
in the ordinary domain where the weight function is defined, that is (0, 00) for the LI and 
L2 polynomials and (—1, 1) for the Jl and J2 polynomials. The behavior of these ordinary 
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zeros are the same as those of other ordinary orthogonal polynomials. In addition to these 
n zeros, there are extra £ zeros outside the ordinary domain. 

For the ordinary zeros, their distribution as n increases at a fixed £ follows the patterns 
of the zeros of the ordinary classical orthogonal polynomials: they are governed by the 
oscillation theorem, and the n + 1 zeros of Pe^+iiv'i ^) interlace with the n zeros of P^ n (j]; A). 
On the other hand, when £ increases at a fixed n, the type LI and L2 Laguerre polynomials 
will have all their n zeros shifted to the right. For the Jl and the J2 Jacobi polynomials, 
both the positive and negative zeros move toward the origin 77 = as £ increases. 

For the £ extra zeros of P^ n {r\\ A), each and everyone of them is sandwiched between the 
corresponding zeros of the deforming polynomials £^(77; A + S) and £^(77; A). As n increases 
at a fixed £, the extra zeros move from the zeros of ^(77; A + 6) to those of £^(77; A). 

The behaviors of the extra zeros as £ increases at a fixed n are more complex. For the LI 
the Laguerre polynomials, all its extra zeros lie on the negative x-axis. So as £ increases by 
one, the number of the extra zeros increases from £ to £ + 1. For the L2 Laguerre and Jl and 
J2 Jacobi polynomials, they have £/2 pairs of complex zeros for even £, and (£ — l)/2 pairs 
of complex zeros and a real zero outside the ordinary domains where the weight functions 
are defined. As £ increases, all the complex zeros move toward the right in the case of the 
L2 Laguerre and Jl Jacobi polynomials, and toward the left for the J2 Jacobi polynomials, 
while the extra real zeros bifurcate into new pairs of complex zeros. For large £, these zeros 
appear to distribute symmetrically with respect to the x-axis in horse-shoe patterns. It 
is interesting to note that in the asymptotic regions of the parameters (g 3> 1, h ^> 1), 
the exceptional polynomial P£ n (i], A) is expressed as the product of the original polynomial 
P n {rf) and the deforming polynomial ^(77; A), (I2"2l . (J25J) and (l2"g|) . 
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Figure 1: LI: Distributions of the zeros %' n \ r\^ ,n ^ (♦), £^ (O) and C) for the 
LI Laguerre polynomials, with g = 0.5 and 1 = 2. The three diagrams correspond to 
n = 1(a), 2(6) and 3(c), respectively. The ordinary zeros rlj'^ lie in (0, oo). 
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Figure 2: L2: Distributions of the zeros fj^' n \ vif'^ (♦), (O) an d (•) for the 

L2 Laguerre polynomials, with g = 0.5 and i = 3. The three diagrams correspond to 
n = 1(a), 2(6) and 5(c), respectively. The ordinary zeros rjf' 71 ^ lie in (0, oo). 
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Figure 3: J2: Distributions of the zeros r]f' n \ rj 

(♦), Ck (O) and Ck (■) f o r the J2 Jacobi 
polynomials, with g = 3, h = 4 and £ = 3. The three diagrams correspond to n — 1(a), 2(6) 
and 5(c), respectively. The ordinary zeros r]f' n ^ lie in (—1, 1). 
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Figure 4: LI: Distributions of the zeros fff' n \ rif'^ (♦), $p (O) an d ^ W f° r the 
LI Laguerre polynomials, with g = 0.5 and n — 2. The three diagrams correspond to 
£ = 1(a), 2(6) and 3(c), respectively. 
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Figure 5: LI: Distributions of the zeros fjf' n \ rjf' 1 ^ (♦), $p (O) an d W f° r the 
LI Laguerre polynomials, with g = 1.5 and n — 2. The three diagrams correspond to 
£ = 1(a), 2(6) and 3(c), respectively. 
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Figure 6: L2: Distributions of the zeros f]f , vf' 7 ^ (♦), (O) an d (•) for the 

L2 Laguerre polynomials, with g = 2 and n = 2. The four diagrams correspond to £ = 
1(a), 2(6), 3(c) and 20(d), respectively. 
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Figure 7: L2: Distributions of the zeros f]f , vf' 7 ^ (♦), (O) an d (•) for the 

L2 Laguerre polynomials, with g = 5 and n = 2. The four diagrams correspond to £ = 
1(a), 2(6), 3(c) and 20(d), respectively. 
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Figure 8: J2: Distributions of the zeros f}^' n \ T]f' n ^ (♦), (O) and (■) for the J 2 
Laguerre polynomials, with g = 3, h = 4 and n = 4. The four diagrams correspond to 
£ = 1(a), 2(6), 3(c) and 20(d), respectively. 
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Figure 9: J2: Distributions of the zeros % , Vj (^)> £fc (O) an d £i W f° r the J2 
Laguerre polynomials, with g = 7, h = 8 and n = 4. The four diagrams correspond to 
£ = 1(a), 2(6), 3(c) and 20(d), respectively. 
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Figure 10: J2: Distributions of the zeros 7] K k 
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